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Abstract 

The paper investigates the synchronization of a network of identical linear state-space 
models under a possibly time-varying and directed interconnection structure. The main re- 
sult is the construction of a dynamic output feedback coupling that achieves synchronization 
if the decoupled systems have no exponentially unstable mode and if the communication 
graph is uniformly connected. The result can be interpreted as a generalization of classical 
consensus algorithms. Stronger conditions are shown to be sufficient - but to some extent, 
also necessary - to ensure synchronization with the diffusive static output coupling often 
considered in the literature. 

1 Introduction 

In these last years, consensus, coordination and synchronization problems have been popular 
subjects in systems and control, motivated by many applications in physics, biology, and engi- 
neering. These problems arise in multi-agent systems with the collective objective of reaching 
agreement about some variables of interest. 

In the consensus literature, the emphasis is on the communication constraints rather than 
on the individual dynamics: the agents exchange information according to a communication 
graph that is not necessarily complete, nor even symmetric or time-invariant, but, in the ab- 
sence of communication, the agreement variables usually have no dynamics. It is the exchange 
of information only that determines the time-evolution of the variables, aiming at asymptotic 
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synchronization to a common value. The convergence of such consensus algorithms has at- 
tracted much attention in the recent years. It only requires a weak form of connectivity for the 
communication graph [H [21 [3l 01 [5] . 

In the synchronization literature, the emphasis is on the individual dynamics rather than 
on the communication limitations: the communication graph is often assumed to be complete 
(or all-to-all), but in the absence of communication, the time-evolution of the systems' variables 
can be oscillatory or even chaotic. The system dynamics can be modified through the infor- 
mation exchange, and, as in the consensus problem, the goal of the interconnection is to reach 
synchronization to a common solution of the individual dynamics [U [H E] . 

Coordination problems encountered in the engineering world can often be rephrased as con- 
sensus or synchronization problems in which both the individual dynamics and the limited 
communication aspects play an important role. Designing interconnection control laws that can 
ensure synchronization of relevant variables is therefore a control problem that has attracted 
quite some attention in the recent years [TQ1 [TJJ [[21 Q31 [J3] . 

The present paper deals with a fairly general solution of the synchronization problem in 
the linear case. Assuming N identical individual agents dynamics each described by the linear 
state-space model (A,B,C), the main result is the construction of a dynamic output feedback 
controller that ensures exponential synchronization to a solution of the linear system x = Ax 
under the following assumptions: (i) A has no exponentially unstable mode, (ii) (A, B) is stabi- 
lizable and (A, C) is detectable, and (iii) the communication graph is uniformly connected. The 
result can be interpreted as a generalization of classical consensus algorithms, studied recently, 
corresponding to the particular case A = [TJ [2|. The generalization includes the non-trivial 
examples of synchronizing harmonic oscillators or chains of integrators. 

The proposed dynamic controller structure proposed in this paper differs from the static dif- 
fusive coupling often considered in the synchronization literature, which requires more stringent 
assumptions on the communication graph. For instance, the results in the recent paper [TS] 
assume a time-invariant topology. The paper also provides sufficient conditions for synchroniza- 
tion by static diffusive coupling and illustrates on simple examples that synchronization may 
fail under diffusive coupling when the stronger assumptions on the communication graph are 
not satisfied. 

The paper is organized as follows. In Section [2] the notation used throughout the paper is 
summarized, some preliminary results are reviewed and the synchronization problem is intro- 
duced and defined. In Section[3]the linear case is studied when state coupling among the systems 
is allowed while in Section 0] the output coupling case is considered. In Section [5] we extend the 
main results to discrete-time linear systems and to periodic time-varying linear systems. Fi- 
nally, in Section [6l two-dimensional examples are reported to illustrate the role of the proposed 
dynamic controller in situations where static diffusive coupling fails to achieve synchronization. 
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2 Preliminaries 



2.1 Notation and Terminology 

Throughout the paper we will use the following notation. Given N vectors X\,X2, ■■■ ,xn we 
indicate with x the stacking of the vectors, i.e. x = [xJ,X2, ■ ■ ■ ,x N ] T . We denote with In the 
diagonal matrix of dimension N x N and we define ljv — [1, 1, • • • , l] r G WL N . Given two matrices 
A and B we denote their Kronecker product with A® B. For notational convenience, we use 
the convention An = In <8> A and A n = A® In- We recall some properties of the Kronecker 
product that will be used throughout the paper 

A <g> B <g> C = A <g> (B <g> C) = (A ® B) <g> C (la) 

A®(B + C) = A®B + A®C (lb) 

iB®CD = (i®C)(B®D) (lc) 

A®B = {A® I p )(I n (g> B) = (I m ®B)(A® I q ) (Id) 

AB®I n = (A®I n ){B® I n ) (le) 

where A G M mn , B G M™. 



2.2 Communication Graphs 

Given a set of interconnected systems the communication topology is encoded through a com- 
munication graph. The convention is that system j receives information from system i iff there 
is a directed link from node j to node i in the communication graph. Let Q(t) = (V, 8(t), Ad(t)) 
be a time- varying weighted digraph (directed graph) where V = {vi, . . . , v n} is the set of nodes, 
£{t) C V x V is the set of edges, and A<j(t) is a weighted adjacency matrix with nonnegative 
elements a,kj(t). In the following we assume that Ad(t) is piece-wise continuous and bounded 
and akj(t) G {0} U [77,7], V/c,j, for some finite scalars < rj < 7 and for all f > 0. Furthermore 
{vk, Vj} G if and only if akj(t) > rj. The set of neighbors of node Vk at time t is denoted by 
Nk{t) — {vj G V : cifcj(i) > 7/}. A path is a sequence of vertices such that for each of its vertices 
Vk the next vertex in the sequence is a neighbor of v^. Assume that there are no self-cycles i.e. 
a-kk(t) = 0, k = 1, 2, . . . , N, and for any t. 

The Laplacian matrix Lit) associated to the graph Q(t) is defined as 



hj(t) 



N 

^2a ki (t), j = k 
-ojfcj(t), j / k. 



The in-degree (respectively out-degree) of node is defined as d\ n = S^=i a fcj (respectively 
d° k ut = J2f=i a jk)- The digraph Q(t) is said to be balanced at time t if the in-degree and the 
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out-degree of each node are equal, that is, 

N N 

^Z a kj = ^2 a jk, k = l,...,N. 

i=i i=i 

Balanced graphs have the particular property that the symmetric part of their Laplacian matrix 
is nonneg ative: L + L T > [16]. We recall some definitions that characterize the concept of 
connectivity for time- varying graphs. 

Definition 1 The digraph Q(t) is connected at time t if there exists a node v k such that all the 
other nodes of the graph are connected to v k via a path that follows the direction of the edges of 
the digraph. 

Definition 2 Consider a graph Q{t). A node v k is said to be connected to node Vj (vj ^ Vi) in 
the interval I = [t a , tb] if there is a path from v k to Vj which respects the orientation of the edges 
for the directed graph (V, Utt=i£(t), JjA^^dr). 

Definition 3 Q(t) is said to be uniformly connected if there exists a time horizon T > and 
an index k such that for all t all the nodes Vj (j ^ k) are connected to node Vk across [t, t + T]. 

2.3 Convergence of consensus algorithms 

Consider N agents exchanging information about their state vector Xk, k = 1,...,N, according 
to a communication graph Q{t). A classical consensus protocol in continuous-time is 

N 

Xk = ^2a kj (t)(xj - x k ), k = l,...,N. (2) 
In discrete-time the analogous dynamics write 

N 

x k (t + l) = x k (t)-e k J2kj(t)xj(t), k = l,...,N (3) 

where e k € (0, l/d k n ). Using the Laplacian matrix, ([2]) and ([3]) can be equivalently expressed as 

x = -L n (t)x. (4) 

and 

x{t + 1) = (l nN - eL n (t)) x(t), (5) 

where e = e ® I n and e = diag(ei, €2, ■ ■ ■ , ctv)- 

Algorithms (Jl|) and ([5]) have been widely studied in the literature and asymptotic conver- 
gence to a consensus value holds under mild assumptions on the communication topology. The 
following theorem summarizes the main results in [T] and [2]. 
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Theorem 1 Let x k , k = 1,2, ... ,N, belong to a finite- dimensional Euclidean space W. Let 
Q{t) be a uniformly connected digraph and L(t) the corresponding Laplacian matrix bounded 
and piecewise continuous in time. Then the equilibrium sets of consensus states of ^> and |5|) 
are uniformly exponentially stable. Furthermore the solutions of fijty and and |5)) asymptotically 
converge to a consensus value ljv ® f3 for some f3 € W . □ 

A general proof for Theorem [1] is based on the property that the convex hull of vectors x k £ W 
is non expanding along the solutions. 

2.4 The Synchronization Problem 

Consider N identical dynamical systems 

&k = f(t,x k ,u k ) (6a) 
Vk = h(x k ), (6b) 

for k = 1,2, . . . , N, where x k € M. n is the state of the system, u k £ M m is the control and 
y k € M p is the output. We assume that the coupling among the systems involves only the output 
differences y k — yj and the controller state differences — According to the graph-theoretic 
description of the communication topology, two systems are coupled at time t if there exists 
an edge connecting them in the associated (time- varying) communication graph Q(t) at time t. 
We will call a control law dynamic if it depends on an internal (controller) state, otherwise it 
is called static. For the systems to be synchronized, the control action (that will depend on the 
coupling) must vanish asymptotically and must force the solutions of the closed-loop systems 
to asymptotically converge to a common solution of the individual systems. This leads to the 
formulation of the following problem: 

Synchronization Problem: Given iV identical systems described by the model ([6]) and a 
communication graph G(t), find a (distributed) control law such that the solutions of (|6|) asymp- 
totically synchronize to a solution of the open-loop system x$ = f(t, xq,0). □ 

In the present paper we focus the attention on synchronization of linear time-invariant sys- 
tems. Generalizations will be the subject of future work. 

3 Synchronization of linear systems with state feedback 

Consider TV identical linear systems, each described by the linear model 

x k =Ax k + Bu k , fe = 1,2, ...,7V, (7) 

where x k € W 1 is the state of the system and u k £ M m is the control vector. For notational 
convenience it is possible to rewrite (|7|) in compact form as 

x = Ajsrx + Bjyu. 
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Theorem Q] can be interpreted as a synchronization result for linear systems with A = and 
B = I. A straightforward generalization is as follows. 

Lemma 1 Consider the linear systems ^Tty. Let B be a n x n nonsingular matrix and assume 
that all the eigenvalues of A belong to the imaginary axis. Assume that the communication graph 
G(t) is uniformly connected and the corresponding Laplacian matrix L(t) piecewise continuous 
and bounded. Then the control law 

N 

Uk = B' is ^a kj (t){x j - xk), k = 1,2, . . . , N, 
exponentially synchronizes all the solutions of ^> to a solution of the system xq = Axq. □ 



Proof: Consider the closed-loop system 

N 



x k = Ax k + '^2a kj (t)(x j - x k ), k = 1,2, ... ,N. 

3=1 



The change of variable 



leads to 



z k — e ^ °^x k: 



k = l,2,...,N, 



N 



N 



z k = -Ae-^-^x, + e-^-^Axk + e ~ A ^ £ a kj (t)( Xj - x k ) = £ a kj (t)(z 3 - z k ) 

3=1 3=1 



or, in a compact form, 



z = -L n {t) z. 



From Theorem [T] the solutions z k (t), k = 1, 2, . . . , N, exponentially converge to a common value 
xq £ R n as t — > oo, that is, there exist constants 5± > and 82 > such that for all to, 



\\z k (t) - so|| < 5^-^-^ \\z k (to) ~ *o\ 
In the original coordinates, this means 



V< > to. 



x k (t) - e A ^-^ 



< <5i 



,A(t-t ) 



- 52{t ~ to) \\x k (t ) - x \\ , 



for every t > to. Because all the eigenvalues of the matrix A lie on the imaginary axis, there 
exists a constant £3 > such that 



x k (t) - e A ^xo < S^-W-rt \\x k (t ) 



for every t > to, which proves that all solutions exponentially synchronize to a solution of the 
open loop system. ■ 
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Remark 1 The result is of course unchanged if A also possesses eigenvalues with a negative 
real part. Exponentially stable modes synchronize to zero, even in the absence of coupling. In 
contrast, the situation of systems with some eigenvalues with a positive real part can be addressed 
is a similar way but it requires that the graph connectivity is sufficiently strong to dominate the 
instability of the system. This is clear from the last part of the proof of Lemma [7] where the 
exponential synchronization in the z coordinates must dominate the divergence of the unstable 
modes of A. 

The assumption of a square (nonsingular) matrix B in Lemma Q] is now weakened to a stabiliz- 
ability assumption on the pair [A, B). For an arbitrary stabilizing feedback matrix K, consider 
the (dynamic) control law 



N 

f} k = (A + BK) r] k + ^2 a kj (*) (Vj ~Vk + x k - xj) 
Uk = Kn k , 

for k = 1,2,... , N, which leads to the closed-loop system 



(8) 



x = Anx + BnKnv (9a) 
77 = (A N + B N K N )r] + L n (t)(x (9b) 



Theorem 2 Consider the system Assume that all the eigenvalues of A belong to the closed 
left-half complex plane. Assume that the pair (A, B) stabilizable and let K be a stabilizing 
matrix such that A + BK is Hurwitz. Assume that the graph is uniformly connected and that the 
Laplacian matrix is piecewise continuous and bounded. Then the solutions of |Sj) exponentially 
synchronize to a solution of the open loop system xq = Axq. □ 



Proof: With the change of variable s k = x k — n k we can rewrite (j9b[) as 

s = A N s - L n (t)s, 

and the closed-loop dynamics write 

x = (An + B n K n Jx + B n K n s (10a) 
s = A N s-L n (t)s. (10b) 

Observe that the two systems (jlOap and (jlObp are decoupled. Since the assumptions of Lemma 
[T]are satisfied for the sub-system (jl0b|) . its solutions exponentially synchronize to a solution of 
so = Asq. The subsystem (f9b|) is therefore an exponentially stable system driven by an input 
L n (t)s(t) that exponentially converges to zero. As a consequence, its solution r](t) exponentially 
converges to zero, which implies that the solutions of (|9ap exponentially synchronize to a solution 
of Xq = Axq. ■ 
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4 Synchronization of linear systems with output feedback 



Consider a group of N identical linear systems described by the linear model 

x k = Ax k + Bu k , (11a) 
y k = Cx k (lib) 

for k = 1, 2, . . . , N, where x k S W 1 is the state of the system, u k S W 71 is the control vector, and 
y k € W is the output. 

The state feedback controller of Theorem[2]is easily extended to an output feedback controller 
if we additionally assume that the pair (A, C) is detectable. Pick an observer matrix H such 
that A + HC is Hurwitz and consider the output feedback controller 

N 

f] k = (A + BK)r] k + '^2a kj (t)(7] j -r] k + x k -Xj) (12a) 

x k = Ax k + Bu k + H(y k - y k ) (12b) 
u k = K7] k (12c) 
y k = Cx k , (12d) 

for k = 1,2, ... ,N, where detectability is assumed and H is a suitable observer matrix. The 
convergence analysis is similar to the one for Theorem [2] and is mainly based on the observation 
that the estimation error is decoupled from the consensus dynamics. 

Theorem 3 Assume that the system ill]) is stabilizable and detectable and that all the eigen- 
values of A belong to the closed left-half complex plane. Assume that the communication graph 
is uniformly connected and the Laplacian matrix is piecewise continuous and bounded. Then for 
any gain matrices K and H such that A + BK and A + HC are Hurwitz, the solutions of Ul\) 
with the dynamic controller Uty) exponentially synchronize to a solution of xq = Axq. 



□ 



Proof: Define s k = x k — r\ k and e k = x k — x k , and rewrite the closed loop system as 



x = \A N + B N K N \x + B N K N (e + s) 

s = Ans — L n s 

e = (A N + H N C N )e. 



This system is the cascade of the closed- loop system analyzed in the proof of Theorem [2] with 
an exponentially stable estimation error dynamics, which proves the result. ■ 
Theorem [3] provides a general synchronization result for linear systems but the solution 
requires a dynamic controller. For the sake of comparison, we provide a set of sufficient conditions 
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to prove synchronization under a simple static output feedback (diffusive) interconnection. These 
sufficient conditions require stronger assumptions on the interconnection and assume a passivity 
property for the system (A, B, C), that is, the existence of a symmetric positive definite matrix 
P > that verifies 

PA + A T P<0, B T P = C. (13) 

Passity conditions have been considered previously in [17] (where it is assumed that the com- 
munication topology is bidirectional and strongly connected) and in [H] (where synchronization 
is studied for a class of (nonlinear) oscillators assuming that the communication topology is 
time-invariant and balanced). Assumptions Al and A2 below lead to a time- varying extension 
of the results in [8] and [17] in the special case of linear systems. 

Theorem 4 Consider system Ul\) with the static output feedback control laws 

N 

Uk = ^2a k j{t){yj - y k ), k = 1, 2, . . . , N. 
j'=i 

Let the graph Laplacian matrix L(t) be piecewise continuous and bounded. Then exponential 
synchronization to a solution of xq = Axq is achieved under either one of the following assump- 
tions: 

Al. The system (A,B,C) is passive and observable, the communication graph is connected and 
balanced at each time; 

A2. The system (A,B,C) is passive and observable, the communication graph is symmetric, 
i.e. the Laplacian matrix can be factorized as L = DD T (t), and the pair (An , Dp (t)C V) is 
uniformly observable. □ 

Proof: Supppose that assumption Al holds and consider the matrix P solution of . 
Consider the Lyapunov function 

V(x) = ^(fl n x) T P N (fl n x), (14) 
the derivative along the solutions of the closed loop system is 

V(x) = -x T fl n P N il n A N x + ]^x T tl n P N tl n A N ±. 



By using the commutation property (jld|) of Kronecker product and the passivity relation ([131) 
we obtain 

v(x) = ^ T n n (p^i 7V + ^^)n„.x-x T ^n p L7 m (t)n p y 

< -y T Ti p Lf m (t)f[ p y. 

Because the graph is balanced, the matrix L sym (t) = (L(t) + L T (t))/2 is positive semi-definite 
for each t and 

(fl p y) T L s p y m (t)n p y > X* 2 tl p y " 
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where X 2 = infj \2(t), and \2(t) is the algebraic connectivity of the graph at time t. Note that 
A2 > because the graph is connected at each time t and the values of the adjacency matrix 
related to the connected components are assumed to be bounded away from zero (see Section 
[2]). This allows to rewrite (|15p as 



V(x) < -X* 2 tl p y , \* 2 >0. 
Integrating (j 16j) over the interval [to, to + T] where T > is arbitrary, we obtain 



(16) 



t +T 



Vdt <-\* 2 



rto+T 




2 








dt < -jX* 2 


tl n x(t ) 


'to 









7>0, 



for all x(to), where the last inequality follows from the observability condition of the pair (^4, C). 
We conclude from a standard Lyapunov argument that the solutions exponentially converge to 
the invariant subspace 



x G R 



nN . 



1 N 



N 



1,2, 



,N 



(17) 



and therefore they exponentially synchronize. To prove that they actually synchronize to a 
solution of the open-loop system it is sufficient to observe that the coupling vanishes in (|17p . 
This implies that the solutions converge to the w-limit sets of the uncoupled system that belong 
to (fT7|) . concluding the first part of the proof. 

Assume that assumption A2 holds. First observe that from the symmetry of the communi- 
cation graph the Laplacian matrix can be factorized as L(t) = DD T (t). Uniform observability 
of the pair (An, Z)JCtv) means that for all to > there exist positive constants T and a (inde- 
pendent from to) such that 



rto+T 
Jto 



^ N (t,t ) T C^D p D^{r)C N ^ N {t,t )dt > aI nN , 

where <I>(i, r) is the transition matrix. This implies that the system 

Amx 



x 

z 



D^{t)C N x, 



(18a) 
(18b) 



is uniformly observable. Applying output injection to system (fT8"|) we obtain 



x 
z 



A N x - K(t)DpC N x 
Dp{t)C N x. 



Choose K(t) = P^Cj^Dp (t) and observe that, since L(t) is bounded, K(t) belongs to L 2 (t,t + 
T). Then output injection preserves observability (see [TB] and references therein) and the 
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system 



x 
z 



Ax - B N DpD p (t)C N x 
DUt)C N x 



is still uniformly observable (here we have also used the passivity condition Cjv = B^Pn)- 
Therefore for all to > there exist positive constants T and (3 (independent from to) such that 
for every x(0) ^ 



f()+T 



\z\\ dt 



to 



t +T 



to 



y(t) T D p f)T(t)y(t)dt > f3. 



Consider the Lyapunov function (|14p . Integrating its time derivative over the interval [to , to + 
T] where T > is arbitrary we obtain 



t +T 



Vdt < 



to 



to+T 2 

Ii p D p Dl{t)y 



dt < -a 



n n x(t } 



a > 0. 



We conclude from standard Lyapunov results that the solutions asymptotically synchronize. 
The rest of the proof is equivalent to the end of the proof under Assumption Al. ■ 



5 Extensions and Generalizations 

In the previous sections the results have been presented for time-invariant linear systems in 
continuous time. For the sake of completeness, we briefly discuss straightforward extensions to 
discrete-time systems and periodic systems. 

5.1 Discrete-Time Linear Systems 

The first step is to provide a discrete-time counterpart of Lemma [TJ Consider the discrete-time 
linear system 

x k (t + 1) = Ax k {t) + Bu k , t = l,2..., k = l,...,N. (19) 

From Theorem Q] we know that the solutions of the system 

N 

z k (t + 1) = z k (t) -e k ^2 hj(t)zj{t), k = l,...,N 

3=1 

where e k € (0, 1/d™), asymptotically converge to a consensus value if the assumptions of Theo- 
rem [1] are satisfied. With the change of variable x k = A^~ to ^z k , t > to, we obtain 

N 

x k (t + 1) = A^-^z^t + 1) = Ax k (t) + ekAj^kji^Xjit), k = l,...,N. (20) 

3=1 
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Identifying (fT9|) and ([20]) results in the control law 

N 

u k = e k B~ 1 AY J h j {t)x j {t), (21) 

3=1 

where we assumed that B is invertible. Assume now that the eigenvalues of A belong to the 
unit circle (in the complex plane). Then there exist 7 > and < q < 1 such that for all to 



< 



soil < 1 q (t ~ to) \\x k {to) - x \\ , k=l,...,N, t>t . 



This shows that the solutions of system (|19|) equipped with (|2ip synchronize to a solution of the 
open-loop system xo(t + 1) = Axo(t). This result is summarized in the following theorem. 

Lemma 2 Consider the system U9\) . Let B be a n x n nonsingular matrix and assume that all 
the eigenvalues of A belong to the boundary of the unitary closed disk (in the complex plane) . As- 
sume that the communication graph Q(t) is uniformly connected and the corresponding Laplacian 
matrix L(t) piecewise continuous and bounded. Then the control law 

N 

u u = e k B- 1 AY,hA t ) x A t ), k = l,2,...,N, e k 6(0,1/4"), 

exponentially synchronizes all the solutions of U9\) to a solution of the system xo(t+l) = Axo(t). 
□ 

Thanks to Lemma [21 we can recast Theorem [2] and Theorem [3] in a discrete-time setting. For 
the sake of compactness we only report the output-feedback case (the state feedback case is just 
a particular case that can be easily derived by the reader). 
Consider the system 

x k (t + l) = Ax k (t) + Bu k (t), (22a) 
y k (t) = Cx k (t), (22b) 



for k = 1, 2, . . . , N, and the discrete-time version of ([12[) 

N 

Vk(t + 1) = {A + BK)r) k {t) + e k Aj2kj{t)(xj(t)-Vj(t)), (23a) 

i=i 

x k {t + l) = Ax k (t) + H (y k (t) - y k (t)) , (23b) 
for k = 1, 2, . . . , N, where y k {t) = Cx k (t). 

Theorem 5 Assume that the system h22\) is stabilizable and detectable and that all the eigen- 
values of A belong to the the closed unitary disk in the complex plane. Assume that the com- 
munication graph Q(t) is uniformly connected and the Laplacian matrix is piecewise continuous 
and bounded. Then for any gain matrices K and H such that A + BK and A + HC are Schur 
matrices, the solutions of IWi) with the dynamic controller [23\) exponentially synchronize to a 
solution of xo(t + 1) = Axo(t). □ 
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The proof of Theorem [5] is straightforward adaptation of the continuous-time counterpart and 
is therefore omitted. 

5.2 Periodic Linear systems 

Periodic linear systems, naturally arise in a number of contexts in engineering, physics, and biol- 
ogy [19] . Periodic models are of large interest also in time-series analysis, economy and finance, 
and in all other cases when seasonal phenomena has to be taken in account. Furthermore they 
arise when linearization of a nonlinear system along a periodic solution is analyzed. Therefore it 
is not difficult to figure out possible applications when synchronization of such models can be of 
interest. The results presented in this section follow from the well-known Floquet theory, where 
the properties of linear periodically time- varying systems are studied via a state transformation 
into a new coordinate system in which the system matrix becomes time- invariant. The eigen- 
values of this matrix are called the characteristic exponents of the original time- varying system 
matrix. In the following the continuous-time case is analyzed, the discrete-time case follows the 
same lines and is omitted for the sake of brevity. Consider the time- varying extension of (|7|) 

x k = A{t)x k + Bu k , k = l,2,...,N, (24) 

where x k £ R n is the state of the system and Uk £ K m is the control vector. The following 
Theorem generalizes Lemma [T] to periodic linear systems. 

Lemma 3 Consider the linear systems {2J$ where A(-) is periodic of period T and let B be a 
nxn nonsingular matrix. Assume that the characteristic exponents of A{-) belong to the closed 
left half complex plane. Assume that the communication graph Q{t) is uniformly connected and 
the corresponding Laplacian matrix L{t) piecewise continuous and bounded. Then the control 
law 

N 

Uk = B~ 1 ^2a kj (t)(x j - x k ), k = l,2,...,N, 
exponentially synchronizes all the solutions of \2J$ to a solution of the system xq = A(£)xq. □ 

Proof: By using Floquet Theory (see for instance [20] ) . there exists a time varying linear trans- 
formation 

so(f) = Q-^xoff) 

where Q(t) is continuous, non-singular and periodic of period T, such that the linear-time varying 
system 

x = A(t)x 
reduces to the linear time-invariant system 

s = ^s , 
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where fiis a constant matrix and its eigenvalues are the characteristic exponents of the original 
system. Moreover the transition matrix can be written as 



$(Mo) = Q(t-t )e 



Cl(t-t ) 



At this point, following the same lines of the proof of Lemma [TJ we observe that with the linear 
transformation = $>(to,t)xk the solutions (in the z coordinates) asymptotically converge to a 
consensus state xq. It follows that there exist constants 6\ > and 82 > such that for all to 



x k {t) - Q(t - t )e n ^x 



Mt-to) 



\\Q(t-to)\\e- 5 ^-^\\x k (t )-x \ 



Vt > to. 



Since Q(t) is continuous and periodic it is also bounded. Moreover the eigenvalues of O (the 
characteristic exponents of A(-) ) are in the close left half complex plane and therefore there 
are no exponentially unstable modes. We conclude that there exist constants £3 > and 64 > 
such that for all to 



Xk{t) - Q(t - t )e 



B(t-t ) 



< 6 3 e- s ^-^ \\x k (t ) 



Vt > t , 



and therefore the solutions exponentially synchronize to a solution of the system xq = A(t)xQ. 



Following the same steps as in Section [31 it is straightforward to translate the results of 
Theorem [2] and Theorem [3] to the periodic case. We leave the extension to the interested reader. 



6 Examples 

The conditions of Theorem|3]are only sufficient conditions for exponential synchronization under 
diffusive coupling. We provide two simple examples to illustrate that these conditions are not far 
from being necessary when considering time- varying and directed graphs and that the internal 
model of the dynamic controller (JSj) plays an important role in such situations. 

Example 1: Synchronization of harmonic oscillators 
Consider a group of N harmonic oscillators 

Xlk = X2k 

X2k = -Xlk + Uk, 

for k = 1,2,... , N, which corresponds to system ([7]) with 

The assumptions of Theorem [2] are satisfied: A is Lyapunov stable and (^4, B) is stabilizable. 
Choosing the stabilizing gain K = (0 — 1), the dynamic control law ([5]) yields the closed-loop 
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Figure 1: The time- varying communication topology used in Example 1 and Example 2. 
system 

X2k = ~X\k ~ V2k 
N 

Vik = V2k + ^ akj(t)(mj - Vlk + xi k - xij) 

N 

mk = -Vik -r}2k + ^2 a fci(*)(^2j ~ ^ 2fc + X2k ~ X2 i)' 

for k = 1,2, ... ,N. Theorem [2] ensures exponential synchronization of the oscillators to a 
solution of the harmonic oscillator if the graph is uniformly connected. Fig. [2] illustrates the 
simulation of a group of 4 oscillators coupled according to the time-varying communication 
topology shown in Fig. [T](the period T is set to 7 sec). The dynamic control ensures exponential 
synchronization. In contrast, synchronization is not observed with the diffusive interconnection 

N 

U k = '^2 a kj(t)(x 2 j ~ X2k)- (25) 

i=i 

The system (A, B, —K) is nevertheless passive, meaning that stronger assumptions on the com- 
munication graph would ensure synchronization with the diffusive coupling f)25[) . We mention 
the recent result [TS] that proves (in discrete-time) synchronization of harmonic oscillators with 
diffusive coupling under the assumption that the graph is time-invariant and connected. The 
following example illustrates an analog scenario with unstable dynamics. 
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Figure 2: First component of the solutions of the closed loop harmonic oscillators by using the 
dynamic control law (to the left) and the static control law (|25[) (to the right). The dynamic 
control ensures exponential synchronization. In contrast, synchronization is not observed with 
the diffusive interconnection. 



Example 2: Consensus for double integrators 
Consider a group of N double integrators 

X\k = x 2k 
X2k = u k , 

for k = 1, 2, . . . , N, which corresponds to system ([7]) with 

'-(::)•■-(:)■ 

The assumptions of Theorem [2] are satisfied: the two eigenvalues of A are zero and (A, B) is 
stabilizable. Choosing the stabilizing gain K = (— 1 — 1), the dynamic control law (JSj) yields 
closeddoop system 

xik = x 2k 

X2k = -Vik - V2k 
N 

rjik = f]2k + ^2 a kj( t )( r llj - Vlk + x lk - xij) 

N 

V2k = -Vlk -1]2k + '^2 a kji t ){ 7 l23 ~ 1]2k + X 2 k ~ X 2 j), 

for k = 1, 2, . . . , N. Theorem [2] ensures exponential synchronization to a solution of the double 
integrator if the graph is uniformly connected. Fig. [3] illustrates the simulation of a group of 
4 double integrators coupled according to the time-varying communication topology shown in 
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Figure 3: First component of the solutions of the closed loop double integrators by using the 
dynamic control law (to the left) and the static control law (|26p (to the right). The dynamic 
control ensures exponential synchronization. In contrast synchronization is not observed with 
the diffusive interconnection. 



Fig. Q](the period T is set to 2 sec). The dynamic control ensures exponential synchronization. 
In contrast, synchronization is not observed with the diffusive interconnection 

N 

U k = '^2 a kj(t)(yj - Uk), Vk = Xik + X 2 k- (26) 

i=i 

The matrix A — aBC is nevertheless stable for every a > 0, suggesting that stronger assumptions 
on the communication graph would ensure synchronization. 



7 Conclusion and future work 

In this paper the problem of synchronizing a network of identical linear systems described by 
the state-space model (A, B, C) under general interconnection topologies has been addressed. A 
dynamic controller ensuring exponential convergence of the solutions to a synchronized solution 
of the decoupled systems is provided assuming that (i) A has no exponentially unstable mode, 
(ii) (A, B) is stabilizable and (A, C) is detectable, and (hi) the communication graph is uniformly 
connected. Stronger conditions are shown to be sufficient (and, to some extent, also necessary) to 
ensure synchronization with the often considered static diffusive output coupling. The extension 
of the proposed technique for synchronization of nonlinear systems is the subject of ongoing work. 
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